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C. Di Troia^ 

^ENEA Unita tecnica Eusione, C.R. Erascati, Via E. Fermi 45, 00044 Frascati (Rome), Ral^ 

The exact analytical description of non relativistic charge motion in general magnetic fields is, 
apparently, a simple problem but, it has not been solved up to now apart for rare cases. The key 
feature of the present derivation is to adopt a non perturbative magnetic field description to find 
new solutions of motion. Among all solutions, two are particularly important: guiding particle and 
gyro-particle solutions. The guiding particle has been characterized to be minimally coupled to the 
magnetic field, the gyro-particle has been defined to be maximally coupled to the magnetic field and, 
also, to move on a closed orbit. The generic charged particle motion is shown to be expressed as the 
sum of such particular solutions. This non perturbative approach corresponds to the description 
of the particle motion in the gyro-center and/or guiding center reference frame obtained at all the 
orders of the modern gyro-center transformation. The Boltzmann equation is analyzed with the 
described exact guiding center coordinates. The obtained gyrokinetic equation is solved for the 
Boltzmann equation at marginal stability conditions. 
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I. INTRODUCTION 


Gyrokinetics is a theoretical framework used for describing either astrophysical or laboratory plasmas behavior 
in the presence of magnetic fields. Its foundation is relied to the motion of a charged particle in a magnetic field. 
The trajectory of a charge is known to be a sort of helix moving circularly around a curve. A point on such curve 
is considered the origin of a spatial reference frame where such gyrating motion is fine described. Gommonly the 
treatment is intrinsically perturbative: the exact solution is known for a uniform and constant magnetic field and 
such solution is thought as the 0*^ order solution. In a seminal paper [I| of the modern guiding center transformation, 
the first equation, here reported for convenience: 


dx^ 

dt 




( 1 ) 


indicates how intrinsic is the perturbative expansion in gyrokinetic theory. The left equivalence is the fundamental 
equation for a dynamical system (without time dependency) and it will be the starting equation also here (with the 
velocity vector field indicated with V instead of A). The right equivalence in (IT|) is always considered in modern 
gyrokinetic theory, expressing a sort of intrinsic expansion useful when the magnetic field is slightly curved and 
stretched and, if time is considered, also time fluctuations are present. The consequence is that particles orbits are 
described with a certain degree of accuracy. Such procedure has been formalized with a huge mathematical effort to 
be extended at the wanted order. The gyrokinetic transformation is such expansion Q. High order expansion respect 
to the leading one must be performed when the magnetic field shows an high level of inhomogeneities in space and 
fluctuations in time. However, this expansion is rarely seen at an order greater than the first. In the present approach 
there is no need of the right equivalence in (HD because the magnetic field, B, is not perturbatively approached. B is 
directly considered nonhomogeneus and fluctuating from the beginning. 

Commonly there is an intermediate transformation introducing gyrokinetics: the Guiding Center (GC) transfor¬ 
mation, which is mainly used when the magnetic field is slowly varying and the time evolution can be neglected (as 
in CD). The reference case is, always, the charge motion in a constant and uniform magnetic field. This example 
clarifies many things and also what the GC is. In a cartesian reference frame the magnetic field is B = BqCz and 
the solution of motion is a velocity decomposed in a parallel and a perpendicular to B component: v = + vj_. 

The perpendicular velocity can be expressed as v± = — |'c±|(ea; siny -|- cosy), where y = Wct is the gyro-angle and 
Wc is the cyelotron frequency, Wc = eB^jm^ where e is the charge, positive for ions, negative for electrons, and m is 
the mass of the considered species (the speed of light is c = 1). The former velocity, v, satisfies the non relativistic 
Lorentz force law: 

V = —{E -I- u X B), 
m 

when, as for example, the electric field E = EqCz is constant and parallel to B. This is easily proved substituting 
■yil = eE^^jm. The motion of the particle is the helix resulting from the uniformly circular motion of the vector Larmor 
radius p around the GG position X which proceeds parallel to B with velocity n||. The particle position is expressed 
hy X = X p, where p = |p|(e 2 ,cosy — Cy siny) has constant magnitude \p\ = pL = ±|i’_l|/wc (the sign depends on 
uJc'- positive if e > 0, negative if e < 0). It is worth noticing that, although it is unstable, the GG itself describes 
a solution of motion and it can represent a real charged particle which is here called Guiding Particle (GP). Such 
particular solution has GG position X and GC velocity V = v\\ez. Note that the GP motion ignores the presence 
of B because it moves parallel to it. If the electric field intensity depends on space, the solution is similar but now 
V = V(x) becomes a velocity vector field. The following three properties will be generalized for several B fields: 1) 
GC can represent a particle obeying the Lorentz force law, 2) thus GP is minimally coupled to the magnetic field, 
and 3) U is a velocity vector field defined at each position (and at all times). 

In the uniform and constant B case, the motion is described in the reference frame centered at the GC position and 
velocity, X and V, respectively. In such reference frame the (unperturbed) trajectory makes a loop, a circular motion 
which is well described by the radial vector p circularly moving with velocity n_L = p. This is again the leading 
example for the gyro-particle motion. The gyro-particle, contrary to the GP, is maximally coupled to B and it moves 
in a closed loop. These two properties will be considered also for the general magnetic field case. 
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For the present studies, it is convenient to slightly modify the formal expression describing the particle velocity. In 
the constant and uniform magnetic field case the velocity is 

V = v\\ez + z;_L = v\\S/z + uic'^S x Vz, (2) 

being Vz = the direction of the magnetic field and, for S = 12, \/S = p, the gyro-radius. The above representation 

can be written as v = U||Vz-|-a;cV x 5'Vz, and considered a sort of Helmoltz representation of v. The first component 
refers to the GP, whilst the last to the gyro-particle, orthogonal either to the GP motion and to VS. Note that 
VS -Vz = 0 then S = v-VS = 0: the surface S' is a constant of motion. The representation © is useful on visualizing 
the helical motion because a helicoid is the sum of two parts: a guiding curve, which for the uniform case is a straight 
line, traced by W|| Vz, and another curve winding around it. If general B is analyzed, the angular velocity vector WcVz 
becomes WcVz — >■ 11 , and the GP velocity v\iVz V, but the same representation allow to recognize the footprint of 
the helicoidal motion. 


II. MOTION OF A CHARGED PARTICLE IN A GENERAL MAGNETIC FIELD 

The most important equation for a dynamical, time dependent, system is 

X = V{t,X). (3) 

If the velocity vector field V is given, the integral curve X could be obtained. Suppose that the field V is unknown 
and you are looking for a solution V which satisfies the non relativistic Lorentz force law: 

V = .^(^E + V X B), (4) 

m 

being B = B{t,X), the external magnetic field, already known in terms of the vector potential A = A{t,X), being 
B = V X A. The electric field, E = —V$ — dtA, being $ = <I>(t, A) the electric potential, is leaved undetermined. 
From the vectorial identity 


V(a • b) = a • Vb -I- b • Va -|-axVxb-l-bxVxa (5) 

applied for a = h = V, the explicit expression of V is 

V{X, t) = dtV + V-VV = dtV + VV‘^/2 -V x{V xV). (6) 

Using the electric and vector potential description of fields in (jl]) and comparing the result with dH]), it is obtained 
the relation: 

dtV + VV‘^/2 -V X {V xV) = - — {dtA + V^ - V x B). ( 7 ) 

m 

Defining the moment per unit mass P = V -\- eA/m, and the energy per unit mass £ = U^/2 -|- e^jm, the above 
equation becomes 


dtP + V£ = V X (v xV + —b) . 

V m / 


( 8 ) 


It follows that 


V -dtP + V -VS = 0. (9) 

The considered system is known to be hamiltonian (the Lorentz^s law is the general Hamilton's equation for a 
conservative system with a velocity dependent potential energy Si) so that it is reasonable setting £ = dt£, or 
V ■ V£ — 0. From dH]), V ■ dtP = 0, or U • dtV = —{elm)V ■ dtA, so that 

£ = V-dtV+ {e/m)dt^ = (e/m) (9*$ - V ■ dtA) . (10) 

Asking for gauge symmetry, P ^ P' = V + eA!jm = V + eA/m + eVgjm and £ —>■ £' = /2 + e^'jm = 

U^/2 -I- e^jm — edtg/m, for a gauge function g. The gauge symmetry ensures the invariance of dH]), being 


dtP + V£ = dtP' -h V£'. 


( 11 ) 
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The equation ([5]) is related to the Euler-Lagrange equation. In the lagrangian formulation, the velocity is an 
unknown variable whilst, in the present approach, it is an unknown function of position and time. There is a 
correspondence between the two approaches due to the fact that the velocity field must satisfies the Lorentz force 
law which can also be obtained from the Euler-Lagrange equation. Analogously, the equation (|S]) is not a Hamilton’s 
equation although the system is hamiltonian. It is worth to note that if you substitute eA/m P —V and e^ jm —>■ 
£ — y^/2 in the Lagrangian, L{X, X,t), for the description of a charged particle in an e.m. field, then L = A^/2 + 
{eAjm) ■ X — e^jm —>■ {X — 14)2/2 -\- P ■ X — £. When the motion is described in the proper reference frame then 
the kinetic energy is zero because the particle is at rest; this case is obtained by posing X = V{t,X), which is the 
equation ([ 3 ]), and the Lagrangian becomes L = P ■ V — £, coherently with the Legendre transformation if the energy 
£ is the Hamiltonian and P is the momentum. 

The equation ([8]) is transformed into the following system of equations: 

idtP-^V£ = -eEt/m, 

\V X (V X V -\- eB/m) = —eEt/m, 


where Et is the transversal electric field. It is an electric field because, from P = V -\- eA/m and £ = 14^/2 + e^jm, 
then 

+ dtV^+Et. (13) 

It is transversal being orthogonal to 14: V ■ Et = 0. Any solution of (USD is solution of motion. The simplest one is 
obtained setting At = 0, or 



dtP + X£ = 0. 

In such case the solution of m satisfies the following fundamental equation: 


^ eB 14 

V X t 4 H-= —. 

TO A 


(14) 


(15) 


The above equation has been already studied in another context, but only if coupled to a second equation for the 
modeling of two magnetofluids although the analogies are interesting, they aren’t considered here. 

The physical meanings of A, which is 


A = 


V2 


{e/m)V ■ B + V -X x 14’ 


(16) 


will be clarified in section|3 In this work, two kinds of solution of equation (11511 are mainly considered, the GP velocity 
and the gyro-particle velocity solutions, corresponding ioV/X'^ (eBjm) and V x 14 ^ —{eB/m), respectively. Before 
discussing the derivation of such kind of solutions, it is worth noticing that equation (fT51) can be rewritten as 


V = XX X P, 


(17) 


which means that a charged particle moves parallel to the curl of the momentum P. In other words, when la is 
verified, a moving charged particle "creates” a magnetic field, B, satisfying m or a vector potential. A, satisfying 
uni). Moreover, from Stake’s theorem, an integral interpretation can be given to (HZD. Suppose to close the orbit of 
the particle in the loop d'E, where E is any normal (smoothed and without holes) surface having 9E as boundary, 
then the flux of 14/A through E is an action. Indeed, 


^v/\ —[ X ' ~ 'f P 'Vdt = J) {£ -\- L)dt 

JT. ^ JaE ./as 


(18) 


being L the Lagrangian. 

In the remainder, the condition (1141) will be applied in (1131) for determining the effective electric field: 


EeS = 


TO 


142 

V— + dtV 


(19) 


When a solution 14 is found for 14 = {e/m){EeS + V x B), then the same solution satisfies 


14 = {e/m) [Aeff + At + 14 X (A + A*)], 


(20) 
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ii Et -\-V X Bt = 0. Choosing Bt = V x Et/V'^, the system (TT^ is rewritten as 


dtP + V£ = —eEt/m, 
y X [V X y + {e/m){B - Bt)] = 0, 


and the solvability equation dia becomes 


V X y + 


e EtxV 

m y^ 


eB 

m 


V 

y 


( 21 ) 


( 22 ) 


As before, A is expressed by (HU), whilst {e/m)Et x VjV'^ is associated, as it will be shown in section 8, to a particular 
’’Ax B” drift velocity. 


III. GUIDING PARTICLE VELOCITY 

The GP solution of (fT5l) is characterized to be minimally coupled to B through the Lorentz force law being V 
mostly parallel to it: V ~ X{e/m)B. This implies that, for a GP, |(ei3/TO)| > |V x y|. The GP velocity is obtained 
substituting: 


V X y = -is • vy, (23) 

where S is a vector of 3 x 3 matrix components. In cartesian coordinates, S is related to the antisymmetric Levi-Civita 
symbol: 




(24) 


When equation m can be inverted, the result is: 


y = (1 + iAS-V) 


i XeB 


being the differential operator (1 + iXS ■ V) ^ defined as 

OO 

(1 + iXS ■ V)"^ = ^(-i)’^A”(S • V)^ 


(25) 


(26) 


n=0 


The equation (1^51) means that the GP velocity belongs from B, after having regularized it with A, and propagating it 
to y. Such solution is the velocity of a particular charged particle, the GP, in a magnetic field B and in the effective 
electric field, consistently given by 


EeS = 


v-^ EdtV]. 


(27) 


It is possible to iteratively obtain the GP solution, (1^51) . starting from the constant and uniform B case. This means 
that also such approach leads to perturbative solutions Q- The leading solution is obtained when B = is 

uniform, and V = is solution of (USD for A = (m/e)v\\/\B\. When B departs from being uniform the velocity 

y(°^ is leaved to be 


y( 0 ) ^ ;^( 0 )£|^ 5 ( 0 )_ 


The first order is obtained inserting the 0*^ solution into the LHS of (|T5]l : 


y(i) = A(i) 


elAl 


6(0) + V X y(°) I = a(i( 

e\B\ 


.e\B\ 


,e|A| 


&(o)+ - vx(A(o)^6W) 

e\B\ m 


(28) 


(29) 


For obtaining A(^) , which is a sort of regularization factor, it is convenient introducing a measurable quantity, namely 
the parallel to B velocity: n|| = y • In such a way, 

A(°) = 


elAl 


(30) 













and 
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, ( 1 ) ^_ {mv\\/e\B\) _ 

1 + {m/e\B\)b^^'> ■ V x (t'||&) 

The first order solution is 

( 1 ) ^ + {m/e\B\)y x {v\\b^°'>) 

1 + {m/e\B\)b^^'> ■ V x (w||&(°)) 


(31) 


(32) 


Such velocity is recognized as the GC velocity (for null magnetic moment) if the GC transformation is truncated at 
first order The velocity at all the orders is simply given by 


ty(") 


+ (m/e|B|)V x 

1 + (TO/e|B|) 6 (o) • V X ■ 


(33) 


Instead of U|| = V-b^^\ it can be used another kind of regularization: the algebraic magnitude |Id|s = sgn(I/ • B)\V\ = 
V ■ b. Note that the unit vector b is parallel to ty, whilst 6 ^"^ is parallel to V^'^\ written as 




6(0) + (m/e|B|)V x ty("-i) 
' '"|6(o) + (m/e|B|)V x y("-i) 


|I"|. 6 (’^). 


(34) 


Gommonly, 6 denotes the unitary vector in the magnetic field direction, but, here, 6 is parallel to V while it is 6 (o) 
which is in the B direction. The advantages of such choice will be clear in section IVIIII and IIXI 


IV. GYRO-PARTICLE MOTION 


In this section it is considered another important solution of equation (1151) which applies when B ^ — (m/e)V x V 
(A —>■ oo). In this case, it looks reasonable taking the vector potential A ^ —{mle)V as solution. This cannot be 
done without conceptual difficulties. The vector potential A has to be defined in the same portion of space where B 
lives and, furthermore, it is gauge symmetric: A ^ A + Vg leaves B unchanged. When A ^ —{mje)V^ both these 
properties seems broken but, only apparently. An example is furnished by the constant and uniform magnetic field 
B = BqBz- For such case, a standard choice of the vector potential is A = rB x Crj^, where the radial vector rCr 
varies in the plane orthogonal to B^ whilst u_l = {e/m)p x B being p = pLCr varies on the circular ring of Larmor 
radius pL- Specifically, A = Air,^) and v± = w_l( 7 ), in poloidal coordinates. The magnetic field B is obtained 
either from i? = V x A as from B = —(m/e)V x v±. Moreover, at r = p^ does exist a gauge function g such that 
= -(jn/e)v^(n) + Vg: 

5=(™/e)^M7, (35) 

being p = {elm)‘^BQp\l2 and V/x = 0; moreover, 7 is the gyro-phase and rVy = e-,. = x Cz- In this gauge, v± is a 
physical, measurable vector potential defined on a ID dimensional space, the ring r = p^. 

For the general magnetic field case, the particle position is a; = A -|- p, and the velocity isv = X-\-p = p = a, being 
X chosen to be zero. In such reference frame, centered at {X,X) = (A, 0), the equation (fTKll becomes 


V X cr -I- 


eB{t, x) 
m 


Kcr, 


(36) 


having replaced V ^ a, X ^ x and A —>■ l/n. The gyro-particle is defined to be the one with V x cr = —{e/m)B(t, x), 
OT K = 0, and whose orbit is closed. As for the uniform and constant case, the gyrophase 7 G 5"^ is the curvilinear 
coordinate describing such closed trajectory. The cyclotron frequency is defined to be the time derivative of 7 : 7 = Wc, 
in this way a = ujcdp/d'j. The ratio |cr|/|a;c| is called Larmor radius which becomes equivalent to pL = Idp/d'yl- 
Applying the Stoke's theorem for computing the magnetic flux 27r'I' of B through the normal (smoothed and without 
holes) surface E, whose boundary, dT, is the closed orbit traced by the gyrating particle, the following result is 
obtained: 


27rT = 


/ B ■ d^x =—{m/e) (j) a ■ dp 
JT, JdT, 




c?7 


= -(m/e) 



(37) 
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or 


/ 2 

(38) 

{eiOc/m) 

So that the magnetic flux through the gyro-particle orbit is expressed by the gyro-average of 27r(m/e)a;cp|, or the 
gyro-average of 2Tr(mle)a^ juJc- There is a particular set of coordinates where such motion is simplified. This happens 
when 7 expresses a local symmetry of B. It will be shown that, in such case, is constant along the trajectory 

and it can be put outside the above integral which becomes: = —(m/e)u}cp\- 

Respect to the general problem of a time varying and non homogeneous it is possible to find a particular 

solution of (ESI) if expressed by 


a = p X V,. (39) 

This solution is similar to the perpendicular velocity in (EJ, with p and 17 to be determined. The vector 17 = 17(7, x) 
is an angular velocity vector field. The curl of a is 

V X (p X 17) = pV • 17 - p • Vf7 -f 17 • V/9 - f7V • p. (40) 

For the uniform magnetic field, 17 is given by 17 = {eBQ/m)ez, and only the fourth term survives. In the non 
homogeneous case it is required that V • 17 = 0 and V • p = I. The relation V • 17 = 0 means that f7’s field lines 
are closed (if not ergodic), as for the magnetic field lines. The relation V • p = 1 means that p describes the spatial 
position of a ID trajectory: p = p(7, 7 ), so that V • p = V 7 • d^p = I. Calling Lpf7 the Lie derivative of 17 respect to 
p: LpLl = p ■ VI7 — 17 • Vp, the curl of a gives rise to the following relation, from (IMl) and (1551) with k = 0: 

eB 

— V X tr = 17 -b LpLl = — . (41) 

m 

If (1411) is inverteble, e.g. for sufficiently small gyro-radius, it could be obtained the solution 

e B 

f7=(I + Lp)-i —, (42) 

m 

being 

00 

(l + Lp)-i = ^(-)"L^. (43) 

n—0 

It is worth noticing that adding to 17 a vector angular velocity along the direction Cp, the gyro-particle velocity in (1551) 
doesn’t change, so that it can be chosen 17 to be perpendicular to p. Another interesting property of such solution is 
the following, being a = p, then 


d p^ 
dt 2 


= p- p = p- a = p- pxfl — 0. 


(44) 


The solution moves on the surface of a sphere whose radius can be identified with the constant Larmor radius pr- 
In this way, the magnetic flux linked to the closed orbit 27r4' becomes, from dSIl), proportional to the period it occurs 
to a particle for one revolution. In other words, magnetic flux can be used as a time-like variable. 

The same solution (1421) can be obtained iteratively, starting from the uniform magnetic field but stretching it a bit 
and leaving it free to fluctuate.. In this way a recursive scheme is found starting from the leading solution: 

= p X 17^°^ = p X {eB/m). (45) 


Applying the curl to the above expression, it gives: 

V X (7^°^ = —eBim 

Multiplying (HHl) per it is possible to define 


- p • V {eB/m) -b {eB/m) ■ Vp. 
as: 

cr(o) . V X cr(°) 


(46) 


( 47 ) 




With the identity matrix 1, and the Lie derivative Lp{eB/m), the equation (j46l) is rewritten as 


V X + eB/m + Lp{eB/m) 



r(0)^(0) 


r(o)|2 




( 48 ) 


The velocity is solution of motion only if Lp{eB/m) ■ = 0. This happens when Lp{eB/m) = 

_fc(o)^(o)^ otherwise next order solution, = p x is processed, can be written as 

— LpLl^^^ = eBim — [p ■ V{eBlm) — {eB/m) ■ Vp]. (49) 


The angular velocity field is divergence free, being Lp{eB/m) 


— {eB/m) — V X The curl of becomes: 


V X = —eB/m + Lp{eB/m) — Lp [{eB/m — Lp{eB/m)] = —eB/m + L^{eB/m) (50) 


or 


V X + eB/m — L^p{eB/m) 





r(l)|2 




(51) 


with, again, 

cr(l) • [V X cr(l) + 

■ 

A recursive scheme is realized: 

{ cr(") = px 

^(n) ^ (_i)"L^(eB/m) 

«:(") = • [V X cr(") + /Icr^^^p. 


(52) 


(53) 


If exists n such that L^f/^\eB/m) ■ [|(T^"^pl — = 0 then is the velocity solution of motion. Otherwise, 

if L^f!/\eB/m) ■ [Icr^^^pl — ^ 0 for all n, and {eB/m)\ < \L/;{eB/m)\, for all n > uq, then the final 

angular velocity vector field, (gH), is the sum of the geometrical series of For such solution, k = 0 and the 

velocity is dMl). 

In the former two sections the time dependency has been kept quite hidden. The problem is that what it has 
been defined to be GP or gyro-particle, can be useful on identifying particles only if it is considered a finite time 
window. Indeed it can happen that for t G At then eB{t)/m ^ V/X, but, before or after a while, for t' ^ At it is 
eB{t')/m ^ —V X V, and V is not anymore a GP. In the following section we describe the generic particles in two 
particular reference frames that take explicitly care of the time dependency, the first is the gyro center reference frame 
(A, V), the other is the GC reference frame that will be indicated by {Q, U). 


V. CHARGE MOTION IN GYRO-CENTER AND IN GC COORDINATES 

The GP and gyro-particle velocities are solutions of the Lorentz’s equation and are obtained from (fT51) . It is 
possible to sum the velocities obtaining a general solution only if the magnetic field is constant and uniform, because 
such solutions live in orthogonal spaces, parallel and perpendicular to B. The aim of this section is to analyze the 
possibility of adding the GP and the gyro-particle solutions to obtain the description of motion for a generic charged 
particle, also for a general magnetic field. The vector position and velocity of the particle, at time t, are written as 

x = X + p, (54) 


and 


V = V + a, (55) 

being]/ = V(t.X) solution of (fTKll . The coordinates (a;, ?;) can be either GP coordinates, if p = 0 and |Vx]/| < |ei?/TO|, 
and/or gyro-particle coordinates, if C = 0 and V x cr ^ —eB/m. As before mentioned, such identification can change 
with time: if |V xV\ < \eB/m\ occurs at time t = 0 but it isn’t anymore true at a later time then the particle is well 
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described as a GP only for a finite time interval near t = 0 G At. 

In this present work the gyro-center is defined by the following statement: the phase space frame of reference where 
the particle orbit makes a closed loop is the gyro-center reference frame Q and indicated by {X,V): gyro center 
position X and gyro center velocity V. It is here assumed that V |t=o= U is solution of m and, moreover, it is a 
GP if p |t=o= 0, with position indicated by X |t=o= Q- The GC coordinates are the coordinates {Q, U), and they are 
defined as the gyro-center coordinates computed at time t = 0 , if it is known that 

= (56) 

m 

being Aq = A(0,(3) and Bq — B{0,Q), the A length and the magnetic field, respectively, computed at t = 0. Such 
magnetic field, Bq, is called equilibrium magnetic field because it is frozen at t = 0: BiBq = 0. The GC coordinates, 
if p |t=o= Oj describe the GP moving in a equilibrium magnetic field. 

The gyro center velocity is not generally known as it is the GC velocity. The difference between X and Q is called the 
GC displacement and it is indicated with The displacement vector depends on the time variation of the magnetic 
field from Bq = B{0, Q) to B(t, X), which is called the magnetic field fluctuation: 5^B = Bit, X) — B{0, Q). It is also 
useful introducing the perturbed velocity, 6^U = V — U, computed by 

C = kU- (57) 

In this way the relation between gyro center and GC coordinates is made explicit: 

X = Q + i, (58) 

V = u + S^u. (59) 

Also the particle coordinates {x, v) can be expressed respect to the GC reference frame: 

X = Q + ^ + p, (60) 

V = U + S^U + a. (61) 

It is worth noticing that only Q and U can be explicitly obtained starting from the equilibrium magnetic field Bq. 
Many kinds of particle motion can be realized with the above description, and it can be helpful consider an orderings 
between such coordinates. It is like for the hands of an analog clock, the hour hand makes a revolution with a slow 
frequency uje, minutes are specified by an intermediate frequency, ujd, and seconds by the highest frequency, lOc, with 
ojc > ujd > ^e- The charge motion is a composite motion in some way similar to the hands motion in an analog clock. 
To specify what time is now, you first look at the hour hand, then, with some approximation, also at the minutes 
hand, and rarely also at the seconds hand. The same is for the charged particle: the GC coordinates can be, like hours, 
the first, rough, representative of the particle coordinates (computed at t= 0 ), the displacement and the perturbed 
velocity can be, like minutes, very useful if summed to the GC coordinates, and a finer motion description is obtained 
by the gyro-center coordinates; the Larmor radius and the a velocity, can be ignored, like seconds. Obviously, this is 
just an example because there are other situations where only the evolution of p and a are considered relevant. 

Concerning the gyro-center coordinates, the particle motion can be easily described starting from the time derivative 
of V, which satisfies the following Lorentz force law: 

e 

V=—{E^s + VxB)=X— + dtV + VxB, (62) 

m I 

being dtP -f Vf = 0 , as seen in (|T^ . 

The energy e of the particle is 

t/2 2 2 

£ = ——I- ^{t, X) H—- —d ■ V -i5p$ = £ -\—-—<j ■ V -\ -(63) 

2 m 2 m 2 m 

being Sp the incremental operator, defined by SpO{t,X + p) = 0{t,X + p) — 0{t,X) for a general field O. If not 
explicitly said, the gyro-center velocity V, the electro-magnetic field E and/or B, and the related potentials $ and 
A, are referred to the time t and position X. The momentum of the particle, p, is the sum oi P = V + [elm)A and 
TT, defined as tt = cr -I- 5pA. Explicitly: 

p = V + a -\ - A H- 6pA = P + tt. 

m m 


( 64 ) 



10 


The Lorentz force law 

V = — [E(t, x) +v X B(t, a;)] 
m 


(65) 


is rewritten as 

a = — [SpE + Vx 6pB] + a x —{B + SpB) 
m m 


( 66 ) 


It is worth noticing that a can be interpreted as the velocity of a particle which is weakly coupled to the ’’electric 
field” x 5pB, that goes to zero together with p. The latter equation has to be compared with the time 

derivative of cr, being a = v{t, x) — V{t, X), so that 


a = dta + vVv-V - VV = dta + V -Va + a -VV + a -Va. (67) 

From the vectorial identity © applied for a = cr and b = (and also for a = b = cr), the Lorentz force law becomes 


dtp + Ve = V X 


V X cr 


-6pB 


+ cr X 


Vxcr + Vxl/ + —B + —SpB 


( 68 ) 


Some of the solutions of the former equation can be chosen to satisfy the following system of equations |ll| : 


f dtp + Ve = 0, 

xV + eB/m = V/X, (69) 

X a -\-V/\ + eSpB/m = 0. 

The velocity cr is ct = p x Lt{t,X + p), as in (IHiHl . even if it it is not anymore a gyro-particle velocity. In this case, 
cr is the particle velocity expressed in the gyro-center reference frame {X, V) at time t. It is worth to note that the 
considered solution is particular being Vxu = VxI/-|-Vxcr = —{e/m)B — {e/m)6pB = —{e/m)B{t,x) 0 - The 
angular velocity n{t,X + p) is obtained from with the substitution eiJ/m —>■ V/X. Explicitly: 

Ll{t,X + p) = {I + Lp) ^ “I-7^)' 

The angular velocity field can be conveniently rewritten as: 

n{t,X + p) =V/X + 6pn. ( 71 ) 

It is possible to give a physical meaning to ujc and A identifying the cyclotron frequency for general B to be the 
algebraic magnitude of Ll(t,x): 


Ll{t,x) = ijJc{t,x)b{t,x). (72) 

If the transit frequency, ujt, is defined to be the limit behavior of Wc for p —0: 

uJt(t,X) = \i-muJcit,x), (73) 

p—tO 

then A can be obtained from: 

= \im^n{t,x) =uJt{t,X)b{t,X). (74) 

The pitch-length A is the path length ran by GP in the time lapse uj~^{t, X). In this way it is re-obtained the known 
property that gyrokinetics is able to describe multi-scale systems. Moreover, if |V x Ej < |ei3/m| the velocity field V 
can be chosen to be 

V ={l+iXS-X)~^—, (75) 

m 

as in (|25|) . In this case the gyro-center velocity V is mostly along the magnetic field direction whilst a is mostly 
orthogonal to B, so that A can be considered the longitudinal scale of the problem, whilst pr the perpendicular 
one. Commonly, both the scales are relevant and the cyclotron frequency determines the ratio u}c{t,x) = \(j\/pr and 
uJcit,x) -)> uJt{t,X) = \V\s/X for PL 0. 
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The problem with gyro-center velocity is that rarely the solution (1751) can be considered correct for a large time 
interval, this happens only for slowly varying magnetic fields. However, the system of equations can be expressed 
by the GC coordinates instead of the gyro center ones: 

( dtp -I- Ve = 0, 

<\7 X 6^U + kU + eS^B/m = (1 - Xok)S^U/Xo, (76) 

I^V X cr -I- (1 — kXQ){U + 6^U)/Xo + e5pB/m = 0, 

where k = 1/X — I/Aq, and U is already known from (1561) . It is worth noticing that the equation for the perturbed 
velocity S^U, as for cr and U, is always the same fundamental equation (fT5|) applied for different ’’magnetic” fields, 
and it is possible to make a group containing all solutions of (1151) . Moreover, it is interesting to notice the hierarchy 
between solutions: U satisfies (USD for B ^ Bq and A — >■ Aq, S^U satisfies the same equation for B —>• {m/e)kU + S^B 
and A —>■ Ao/(l — kXo), or, in other words, the solution for S^U is slave of the solution for U and Aq. Once obtained 
the gyro center velocity V = U + 6^U, the same remark is true for cr, which is solution of (ITSl) for B —>■ {mje)VjX, or, 
cr is slave of V and A solutions which, in turn, are slave of U and Aq solutions, which in turn, depends on Bq. 

In summary, the considered particle is represented in gyro-center coordinates by 

x = X + p, (77) 

v = V + pxil = V + pxV/X + px i5pH, (78) 

where {X, V) are computed at each time t in the gyro-center description or substituted to X = Q-\-^ and V = U + S^U 
in the GC description. 

Up to know it has been considered general magnetic fields and two interesting solutions for the particle motion, in 
the next section it is considered a particular magnetic field, often encountered in the description of laboratory and 
astrophysical plasmas, and the corresponding general solution for the particle motion. 


VI. CHARGE MOTION IN AN AXISYMMETRIC MAGNETIC FIELD 

In this section the charged particle motion is considered when the magnetic field shows some symmetries. If the 
system is symmetric in the toroidal (j) direction, the magnetic field can be written as 

H = VV' X V((i -b FV(/>, (79) 

being (j) the toroidal angle, = e^/R (e^ the unitary vector in the 4> direction and R the major radius in cylindrical 
coordinates), tjj the poloidal magnetic flux. Moreover, V'0 • X/cj) = 0 and d^F = 0, for the considered symmetry. This 
representation, mostly used for describing tokamak magnetic fields, is very similar to the one used on describing the 
velocity in equation ([2]). Indeed, it represents a magnetic field made by helicoidal magnetic field lines and can be 
decomposed into the sum of a toroidal field, Bt, whose direction is along the magnetic axis, and a poloidal field. 
Bp = ±\Bp\ee, for a poloidal angle 9 in toroidal coordinates. Magnetic field lines can be thought as proceeding towards 
the magnetic axis and simultaneously winding around it. Such helicoidal magnetic field lines are characterized by 
the winding number is = q~^, being q the safety factor. 

For a closed magnetic field line, n toroidal cycles match m poloidal ones, so that q = m/n is & rational number, 
q G Q. There is a flux coordinates system where such description is easily visualized because the magnetic lines are 
mapped into straight lines. The surface made of magnetic field lines corresponding to a given q is identified by the 
value of tf, proportional to the poloidal magnetic flux (j)p = 2Tr'ip. If q G Q then ij) is said to be resonant. It is useful 
introducing a (toroidal) radial coordinate which labels tf, in this way r = r{^p) and q = q{r). The radius r depends 
on the intensity of magnetic field, for this reason, r is considered a generalized radial coordinate. It is also introduced 
a generalized poloidal angle, d G S^, so that: 

1 _ B-Vd 

q{r) B-Vf 

In the straight flux coordinates the magnetic field is expressed by 

B = Xfjx X{(j)-qd). 

From dSID, the vector potential is immediately found: 


(80) 

(81) 


A = t/jX/cj) + q-dX/fj + X/g. 


( 82 ) 
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Thus vector potential representation is multivalued. However, it can be shown that 2TTq'V'ijj is the gradient of the 
toroidal magnetic flux (pt = —27rtE' that becomes the gauge function that leaves B single valued. Moreover, such 
representation shows that ip can be defined through the toroidal component oi A: ip = |V(^|“^(H — Vg) ■ e^. A 
further simplification on fields mapping is obtained with the following Clebsh potentials: 7 and'!'. They are defined 
as 7 = — (p/q and from V'l' = —qVip. The magnetic field can then be written in the Clebsh representation as 

B = V'^ X V 7 , (83) 


while the vector potential becomes: 


A = 4'V7 + V(7, (84) 

and can be defined as = |V 7 |“^(A — V 5 ) ■ e^. Note that the Clebsh representation can be adopted also when B 
has no symmetries p^ . 

Armed with the above magnetic field representations it is now possible to address and, in some circumstances, to 
resolve the general case of axisymmetric system. Firstly, equation (EH), for representing B, is used into the equation 
of motion m- The velocity v, solution of m is decomposed in the sum v = Vp + vq. The particular solution, Up, 
satisfies 


eB Vn 

V X 7;p H-= —, (85) 

m Ap 

while Vo is the homogeneous solution satisfying: 

V X no = ^. (86) 

Ap 

Note that (I5S1) is known as the curl spectrum equation, or as the force-free equation [3) Hi • uo is often called a 
Beltrami field p^ . It has been used in p^ to find a magnetic field with fixed helicity and minimum energy for the 
description of some kind of astrophysical sources. As for example, it is reported below the solution find by [13 for 
the spheromak, which is the general solution of equation (I86L when A is constant, in spherical boundary domain, 
Cr -vq = 0, for Cr orthogonal to the ball of radius a. Writing vq = UrCr + ugee + u^e,/, in spherical coordinates (r, 0, <p), 
then 


= {l^-\)r-hi{r/X)P:Picos9) 


COS m(t) 

sin m- 


(87) 


ue = -dr[VrJi{r/X)]dePPI^ {cos 9) 

r 


COS mcj) 
smmd 


T 


m 

\x/r 


jVpM) Cr(cos 6 >) /sinmcp 

oin \ , 


nr\c m / 


( 88 ) 


77 ? 

Ucp = T—dr[VrJi{r/X)] — 


^{cos9) /sin m(p\ m 


sin 9 \ cos md) 


- Y^7K’’/^)'9eCr(cos6») 


cos m(p 


(89) 


with / = n + 1/2, m (not to be confused with the mass) is m G {0,1,..., n} and Ji are Bessel functions (of first kind, 
see (11461) 1 and P™ are associated Legendre functions. 

Coming back to the toroidal symmetric system, Vp and Xp can be easily chosen to be: 


Vp = —{e/m)ipVcp and Xp = —ip/F. (90) 

Such particular solution is proportional to the toroidal component of A\ Vp = —{e/mfA^eci,, see (|H^ with d^g = 0. The 
Beltrami field vg is obtained starting form the following representation (to be compared with the velocity representation 
in (El): 


uo = X Vcp-h PoVpj, (91) 

where the quantity in parenthesis is a special magnetic field with the poloidal magnetic flux 2ttP^ and a toroidal 
component Hg/R, being |V^| = R~^. Moreover, for preserving the toroidal symmetry d^pPcj, = drfjHg = d^vg = 0. 
The curl of (IMl) gives: 


V X Uo = Vt'o X (VP^ X \7(p) + l/gV X W(j, X \7(p + S/{ugHg) X V(/. 


(92) 
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Being 


Vi'o X (V7^0 X V(/)) = -(Vt'o ■ VP^)V4> 


and 

V X (VP^ X V<^) = -p2V(^ V • ^ = -V<^A*P^, 

where A* is the Shafranov differential operator. The equation m becomes: 

Fu 

- (Vi/o ■ VP0)V</) - V<^A*P^ + V(!/oJ^o) X V<^ = —^[VP0 X V<^ + PoV<^]. (93) 

Comparing the toroidal components: 

V^o • VP^ + A*P^ = (94) 

and the poloidal ones: 

Fiy 

V{voHo) xV4> = —^VP0 X V4>, (95) 

The following solution for vq and Pq, are obtained: 

Eih FFa^ 

vq = = {e/m)Xp and Hq = - — = V^l\. (96) 

The total velocity is 

V = Vp + VQ = - X V()| H-(P^ - 'ip) V(p- (97) 

mh m 

The toroidal component of v clarifies the meaning of P^, being proportional to the toroidal momentum = Cz-xxp = 
Cz ■ X X [u + (ejm) x A]: 


Vfj,='ip + —Rv^ = mp^/e, 


(98) 


In an axisymmetric (and isolated) system is a constant of motion, which is immediately verifyed for P^, being 


n • VP^ = 0, (99) 

directly from the chosen representation dUD of vq. The velocity is expressed through a special magnetic field (muJt/e)b. 
Indeed, v = (raje)uit\ph, coherently with equation (|7i)) . being h parallel to v and (vn/e)uith = VP^ x V(/) + PV(^, 
with P = P + V^/Xp = P(I — V(j,/ip). As done for the magnetic field, the velocity field lines lie on the magnetic flux 
surface P^ and a velocity wmdmg ’number, i^, is defined to be 


_ u-Vd 

"" vX/(p' 


( 100 ) 


analogously to the inverse of the safety factor. When i^ S Q then P^ is called resonant toroidal momentum and 
particles trace closed trajectories. It is worth noticing that for closed orbit the linked flux of v/Xp is related to the 
integral in (fT^ . Moreover, if the orbit is closed the frequency to make a cycle in (p is commensurable with the 
frequency uib to make a cycle in d, and both are multiple of Wj-es, the characteristic frequency, or resonant frequency, 
of the closed orbit. Moreover, for the closed orbit, = Wb/wd G Q and WresjWd and Wb are linearly dependent over Q. 
Analogously to the resonant ip, the resonant P^ plays a fundamental role in the gyro-kinetic dynamics as it can be 
seen in the analysis of Poincare maps studied in [18[. 

The toroidal momentum in flux unit, can be obtained, once F is known, from the explicit solution of equation 
(IMl) . here concisely rewritten: 




XpR^ 




(101) 
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being Ap = —ip/F with the operator As = SV ■ equivalent to A* if S' oc R^, which means that Ap is a constant, 

OT F (X Ip. 

Two solutions for without solving equation (I101|) . are described here, for the GP velocity V, and in the next 
section, for the gyro-particle velocity, a. From equation (1551) . the GP toroidal momentum is 


in) ^ I 1 + imle)iRlF)e^ ■ V x mFu|| 

^ 1-h (TO/e)( 6 (o)/|B|) • V X e\B\ ’ 

so that = Ip + {m/e)Fv\\/\B\. Otherwise, from equation (IMl) . becomes 

(„) ^ 1 + im/e){R/\B\)e^ ■ V x mF\V\s 

^ |l-h(m/e)(5(o)/|S|)-VxFO-i)| e\B\ 


( 102 ) 


(103) 


with the new regularization schema proposed (|F|s instead of U||). For the GP velocity V ^ —{e/mPipB/F, the 
velocity (l97ll can be approximated to 

- X v</. + - (P^ - V') yep -^y^ xVeP- (104) 

mr m mr m 

which means that ^ 0 and, yP^ ^ yip or F —>■ oo. Indeed, the GP velocity is mainly directed towards B. The 
drift velocity expresses the deviation of the GP from X(e/m)B. In the present work, the drift velocity is indicated 
with Vd and defined as: 


VD = \y xv = --^y{VA-ip)xV(p + ^V(P, (io5) 

toF m 

where the last equivalence is obtained for A = Ap = —ip/F and V = {e/m)XpS/V^ x \7(p + {e/m){V^ — ip)ycp. It is 
worth noticing that wy = V x F is sometimes called the vorticity, so that Vd = XpLOy- 

It can be consistently checked that Vd is smaller than V when ^ 0 and, ~ yip or F —oo. It is worth 

noticing that the last condition means that the magnetic field is mostly toroidal, as it happens in tokamaks where 
GPs are believed to have an important role. 

From the definition of drift velocity, it is also possible to obtain another useful expression of F^, multiplying both 
sides of (11051) per (m/e)Fe 0 we obtain: 


^ y (106) 

e eF 

The last equation, together with (IIOII) allows to express the gradient of for the general case (F —>■ v, which is not 
anymore the GP velocity): 


VF0 


mXpRf 

e 


V X Vcp. 


(107) 


VII. GYRO-SYMMETRY 

The conserved momentum for the gyro-particle velocity cr, when X = 0 and V x cr = —{e/m)B, needs another 
treatment. The representation (IHIT)) is better conceived for the present case because it describes general magnetic 
configurations and, specifically, also not symmetric magnetic fields. The particular solution of (I36|) . like in (19011 . is 

CTp = —(e/m)'I'V 7 for Ap —>■ oo ( or /Cp = 0). (108) 

The homogeneous equation is reduced to V x o-q = 0 with the solution (Tq = —(e/m)Vg so that 

a =—{e/m){A + y g), (109) 

This result was already shown for the uniform magnetic field case. As in (1551) . a good gauge choice is the following: 

g = {mlef)i'y, (110) 

if g is constant: V/r = 0. The gyro-particle velocity becomes 


a = —(mje) [/r -I- (e/TO)^4'] Vy. 


(Ill) 
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The same expression, (llll|) . just obtained without requiring any symmetry, can be also deduced from (j97l) . when a 
symmetry is evident. With the following substitution: (?!)—> 7 , '0 and —>■ —(m/e)^/r, the velocity is 

a = -^ ^7-+ (e/mV'^] V 7 ; (112) 

er e 

for the gyro-particle Ap —>■ 00 or F — 0 , so that Vfi = 0 must be imposed and the relation (lllll) re-obtained but, now, 
7 is an ignorable variable expressing the symmetry of the system and fj, is the constant of motion associated to such 
symmetry. This means that the gyro-symmetry in 7 was also present, but hidden, in the first derivation of (dm). 

In analogy to the magnetic moment is defined as (m/ej/i = —b ■ p x [a (e/m)A], for the unit vector b = fl/coc, 
where the angular velocity vector is defined in and, for the symmetric case, d-^ojc = 0. This definition of the 
magnetic moment for the gyro-particle is consistent with dm if V 7 • 6 X p = 1, or V 7 = e-y/pi and e-y ,6 and Cp 
define a tern of orthogonal unit vectors basis. 

The gyro-symmetry is important because is linked to the property of the gyro-particle to be closed, that is by definition 
in the present approach. In other words, when there is a closed orbit, it is possible to define a magnetic moment, a 
gyro-angle and a gauge, product of the two times (m/e)^, so that such particle feels the force on itself independently 
from the gyro-angle. In this way it is assumed possible to have Wc = u!c(t,X, p^) or 

= 0 (113) 

also for general magnetic fields. This assumption is relevant because, from ([m, dc = ^{t,X,pL) = -{rn/e)uicp\- 
Moreover, from the representation (1831) . T can be identified with the toroidal magnetic flux, and the Larmor radius 
Pl, can be identified with the radial coordinate used to label the magnetic flux surfaces, pl = f (see the introduction 
of section EH- If X lies on the magnetic axis, it seems reasonable identifying the gyro-angle j to he j = -d — (j)/q, as 
if considered a generalized poloidal angle for the magnetic field expressed in toroidal flux coordinates. 


VIII. MAGNETIC MOMENT AND DRIFT VELOCITY 


In electrodynamics the magnetic moment p, is considered an adiabatic invariant of motion. Such crucial concept 
is in conflict with the present purpose of obtaining a non perturbative description of motion. Indeed, the adjective 
’’adiabatic” refers to an intrinsically perturbative procedure, again. The magnetic moment for a gyro-particle has 
been already considered in (dm, explicitly rewritten below 

p = —b - p X [cr -I- {e/m)A] = —{e/m)p\(T ■ V 7 — (e/m)^d>, (114) 


as the conserved momentum associated to the gyro-symmetry, which allows the orbit of a gyro-particle to close. For 
the charged particle expressed by dm, p can be similarly defined as proportional to the conjugate momentum (in 
hamiltonian sense) of the ignorable coordinate 7 . 

The cyclotron frequency has already been defined as Wc = 7 , with d-fiOc = 0. 

The Hamiltonian, H, associated to the description of the particle motion is 


, b - (e/.n)(,4 + S,A)r ^ 

z m 

= ^ + -+-$ + Sp^) 

2 2 \ m J m 

It is now possible to define the magnetic moment as 

ecr^ ea ■ V / e \ ^ (5„d) 

M = 7^ -^-^ (“) -’ 

ZmcUc mujc \m/ lOc 

so that the Hamiltonian becomes 


%{t,x,p) = H{t,X,P, (m/e) 7 ,/r) = 


[P — {e/m)A]'^ e 


m 


— $ H-^cjc- 

m e 


(115) 

(116) 


(117) 


(118) 


and the motion is described by the canonical coordinates {X,P) and {{m/e)j,p). The Hamiltonian equations for the 
last couple of coordinates are trivial: 


= dpH = {mle)uic, 
p = —{m/e)djH = 0, 


(119) 
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with the desired conservation rule. Moreover, for the coordinates {X,P): 

{ X = VpH = P-{elm)A, 

P = -VH = {e/m)[P - {e/m)A] ■ VA + {e/m)[P - {e/m)A] x V x y4+ (120) 

— (e/m)V$ — {mje)^yu}c- 

The last equation, after some algebra, is nothing else that the Lorentz force law with an effective electric field given 
by Ees = E - (m/e)^Va;c: 


V = -{E,s + V X B). (121) 

m 

There is not any canonical neither non canonical transformation from {x,p) —>■ {X,P,"f,fj,) because the dimension of 
the two spaces is different. Fortunately, this can be done without problems advancing P to be a held P = P{t, X). In 
this way, also V becomes a held V = V{t,X) and the hrst Hamilton's equation in (11201) gives the starting equation 
dSl) of the present work. 

When V ^ eB/m, the drift velocity, Vd = AV x V, if explicitly computed, shows we are in the right direction. 
Starting from the curvature vector Kb of F, being Kb = h ■ V& = —b x V x 6, then 

V X b — b X Kb + bb -V x b. (122) 


The drift velocity becomes: 


Vd = X {V\V\s xb+ |F|s {bx Kb + bb-X xb)} . (123) 

which can be rewritten, for F/A = ojtb, or |F|s = Awt, and F^/2 = e — (e/m)<l> — {m/e)pLUJc, as: 

Vp = X rv$ + nX(muJc/e) + V^Kb] + —V -X xb. (124) 

UJt UJt UJt 

Substituting Ve = —dtp = —dt[p— {e/m)A] — {efm)dtA, Vd becomes: 

Vd = 7 — -r^ -h —& X X{muJc/e) H- b x Kb -\ - b x dAp — {elm')A\ H-F - X x b. (125) 

(mle)LOt oJt oJt oJt oJt 

The drift velocity Vd is the sum of hve terms, respectively called: the E x P-like drift, the /xV|P|-hke drift, the 
curvature-like drift, the 9t|P|-hke drift and the Banos term. 

The velocity Vd is identical in form to the drift velocity computed at the lowest order in standard gyrokinetics 0 - 
Respect to the standard gyro-center drift velocity, e|P|/m must be replaced by Wc and/or by Wt, whilst the unit 
vector in the B direction, must be substituted with b, parallel to F. This means that all results obtained up to now 
with a drift velocity computed at lowest order, can be easily extended and considered almost right at all orders if you 
give a different interpretation to the ambient magnetic field, as if you are considering another scenario, surely not so 
different, but another one. Such remark refers mainly to the laboratory plasma modeling, where there is an optimal 
control on the equilibrium magnetic field; e.g. in fusion devices, where, among other, the configuration of the ambient 
magnetic field is fundamental for the stability properties of plasma and its confinement. 

It is worth noticing that it is possible to obtain the SjlPj-like drift considering Ve = 0 and adding to the drift 
velocity the term b x dtpjiot- This procedure is the same outlined in section HIl when it is considered an electric 
field given by P = Peff + Et and, for the present case, Et = {e/m)dtp. Setting Ve = 0 means considering X and e 
as independent coordinates. There is an advantage on describing the system through the coordinates {t,X,e, p,j). 
Indeed, the charge in a general magnetic field can be described either in canonical coordinates {x,p) as in non canonical 
coordinates (X, e,/x, 7 ). However, those coordinates can be, sometime, recast into canonical one as in the toroidal 
symmetry case. In such case, it is possible write X as function of X = X{tp, P^, {elm)(f>), as done in and the motion 
can be described by the canonical pairs of variables: {t, e), ((e/m)(/, P^) and ((m/e) 7 , p); considering ^ as a time-like 
coordinate, as precedently suggested, and the conjugate to momentum, p^ = p^{t,{e/m)(j),(rn/e)^,£,'P^, p), as 
Hamilton^ function. Such canonical description is postponed to a future work. However, the same set of coordinates, 
(t, X, e, p, 7 ), is used for obtaining the gyrokinetic equation, in the next section. To proceed and obtain the gyrokinetic 
equation it could be helpful to explicitly express the magnetic moment and the drift velocity for a constant magnetic 
field, in GC coordinates: 


Mo 


ecr^ ea ■ U ( ^ \ ^ 

2mu!co rnwco uico ’ 


(126) 
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and 


respectively. 


Ud 


b X V<i>o /To , , N , 17 ^ , 

7 —-1- b X VlmuJco e) H-5 x H- U • V x 6, 

[m e)wto wto Wto Wto 


(127) 


IX. GYROKINETIC EQUATION 


The gyrokinetic equation is chosen as the first use of the present approach. The gyrokinetic equation is a fundamental 
equation for describing plasma behavior in a magnetic field. Up to now, it has been studied the single particle motion 
in general magnetic fields, whilst the gyrokinetic equation analyzes a macroscopic system and it is inherent to the 
many-body physics. Thanks to the Boltzmann equation you can distinguish what depends on the dynamic of the 
single charge from the rest, The distribution function and the collision operator, C. The collision operator applied to 
the distribution function / = f(t,x,v) can be chosen to be the modified Landau collision operator: 


Cl = 




1 - 


{f{t,x,v')yyf{t,x,v) - f{t,x,v)Vy>f{t,x,v')). 


(128) 


being yc a constant which is known for the Coulomb collision case, and Up given by 


Up = V — v' + 


pT 

OJ- on 


mA?, 


a 2 {pfi - s§) - A2 cr' 2 (/r'2 - sg) - A^ 


Wc M + So w' n' + So 

to allow the system relaxing towards the following equilibrium distribution function | 


feqiV^, e, p) = N exp 


r in-no?] 


e 



<1 

1 

_1 

exp 

7^ 

w . 

exp 

[ J 


(129) 


(130) 


being A, sq, A^, Tu,, and Ap^ constants parameters. 

However, hereinafter, C and feq are left undetermined. This means that the gyrokinetic equation we get to show 
cannot be exactly solved, here. The main difference with other derivations is that the gyrokinetic equation is here 
obtained non-perturbatively, which means that the single particle motion is described without approximation, using 
the exact magnetic field. Perhaps inevitably, the solution of the gyrokinetic equation have to be expressed perturba- 
tively. 


The form of gyrokinetic equation depends on the set of coordinates used to represent the charge motion. The 
gyrocenter coordinates t, X, e, p, and 7 are the present choice. The time derivative of those variables are already 
known, from ([3]), (fTUl) and (II19^ : 

X = V, £ =— [dt^{t,X + p) - iy + a) ■ dtA{t,X + p)], /i = 0, i=ujc- (131) 

m 

The Boltzmann equation for a distribution function / is 

f = c{f), (132) 

or, more explicitly, 

dtf + V-Vf + - [dMt, X + p)-iV + a)- dtA{t, X + p)] dj + = C{f). (133) 

m 

The initial distribution function is assumed to satisfy Boltzmann equation: /o = Clfo). For stability studies, it is often 
considered an equilibrium distribution function, like feq in (I130L which is a particular solution of the Boltzamann 
equation, being feq = C{feq) = 0, separately. The evolution of plasma is obtained by the behavior of the distribution 
function and the behavior of fields. The Boltzmann equation gives rise to an initial value problem: it is studied the 
system evolution from an initial state, the one defined at time t = 0. All the increments, of the distribution function 
and fields, are referred to the initial values as below explicitly written: 

f{t, X, v) = /o(0, Q, £ 0 , po) + 6f{t, X, £, p, 7 ), 

$(t, x) = $ 0 ( 0 , x) + d$(t, x), 

A{t, x) = Ao(0, x) + SA{t, x), 


(134) 

(135) 

(136) 
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being Q the GC position, and Eq = C/^/2 + e^ojm + {m/e)nQU]cO, the energy at t = 0, with U the GC velocity, /xq 
and Wco, the magnetic moment and the cyclotron frequency, respectively, always at t = 0. It is here assumed that 
the initial distribution function satisfies the Boltzmann equation: /o = C{fo) and that d-yfo = 0. Due to the latter 
condition, /o is said to be the initial distribution function of GGs (instead of particles). The Boltzmann equation for 
/o is particularly simple: 


t/-V/o = C(/o), 


(137) 


being dtfo = do = Mo = 0. 

The gyrokinetic equation describes the functional form of Sf in terms of initial fields and the initial distribution 
function, supposed given. Using (11371) . (11331) and (11371) . the Boltzmann equation is rewritten as 

dtSf + U • Vd/ + LOcd^Sf + S^U ■ V/o + -idt6<S> -V-OtSA-a- dtSA)dJo = 

m 

= 6sfC- -{dtS^ -V-dt6A-a- dtSA)deSf, (138) 

m 

being SsfC = C{f) — C(/o). In this way it has been moved on the RHS the non-linear terms leaving on the LHS all 
the linear terms. The RHS can be computed only if C is substituted with an analytical operator, as can be Cl in 
mB, but here C is leaved undetermined and the RHS is indicated as " —vneST' , being a symbol. Explicitly: 

dtSf + V ■ V5f + uJcdjSf + S^U ■ V/o + ^ ‘ ^tSA - a ■ dt5A)defo = -vnl^!- (139) 

The equation p39D is the Boltzamann equation for Sf which represents all the charges in the distribution function 
that, at time t 7 ^ 0, cannot be anymore represented by /q. Such equation is studied in the dual Fourier space. The 
Fourier transform in (t, X, 7 ) of 5f is Sfi defined as: 

00 p 

Sf^Yl / (140) 

I —— 00 

with the function fi{uj,k) depending also on energy, magnetic momentum and GG displacement vector, even if not 
explicitly indicated. The Fourier transform of e.m. potentials are: 


5<^{t,X + p) = 

J d^kduje^‘^*e-^'^-^^+p'^S^{u;,k), 

(141) 

5A{t^ X -\- p) = 

J d^kduje^‘^*e-^'^<^+P^SA{u;,k). 

(142) 

At least, the Fourier transform of S^U is 




d^kdu;e^‘^*^e-^^-^S^U{uj,k). 

(143) 


The equation (11391) . in the Fourier dual space, is 

dt5f FV-^/Sf + LUcd^Sf + VNL^f + 

+5^U ■ V/o -I- — - V ■ dtSA - a ■ dtSA)dsfo = 

m 

00 p 

= iYl / WArL)d/i-h 

I —— 00 

+ J d^kdoje^^*e-^'^-^S^U ■yfo + 

+i [ d^kdu;e^‘^*^e-^’^-^{e/m)u;e-^^-P{6^-V-SA-a-SA)d,fo = 0. ( 144 ) 


The same equation is below written with the substitution a = uJcP x b, being djUic = 0, from (firm . and the tern of 
orthogonal unit vectors Cp ■ b x e^ = 1 (computed at time t and position x = X + p). When equation (11441) is anti 
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transformed in 7 , then 


/ + V-V6f + ojcd^Sf + -{dtS^ -V-dt6A-a- dtSA)d,fo - i^nlSJ] = 

J Z'K m 


d^kduje^^^e 


iujt ^ — ik‘X 


{uj-k-V - I'wc - ivNL)Sfv +SroJ dPkduje^'^*e-^'^-^6fi ■ V/o + 


+i / d^kduje^^^^e 


iojt^-ik-X 


euj 


[(<5$ -V-SA)- uj,b X SA ■ Vkjdefo 


^^ilA-ik-p _ Q 
2TT 


(145) 


The wave number k varies in the dual space but it can be conveniently chosen to be expressed in cylindrical coordinates: 
k = ki,b + kt = kijb + |fct|[ep sin (7 — 'jk) + e^ cos (7 — 7 ^)], being b the longitudinal (parallel to 17) unit vector, fcf, the 
longitudinal wave number component, \kt\ the amplitude of the transversal wave number component and 7 ^ an 
arbitrary gyro-angle. The cylindrical Bessel functions, Ji, are defined by 


Jii\kt\pL) = ^ f '^9'® 

whose partial derivatives respect to k are 


1 P p — illk 

= _ p^pil{l-7k)-i\kt\pLsin{'y-jk) _ 7 _ 

27r 


d'ye 


il-y—ik-p 


ykM\kt\PL) = yk\kt\■ di,,iM\kt\pL) = ^Jli\kt\PL), 


and satisfying the following property: 


Ji = 


(Ti_i — Ji+i)/2, if 1 ^ 0 
-Ji, if 1 = 0. 


(146) 


(147) 


(148) 


(149) 


Equation (114511 . with the use of the cylindrical Bessel functions, is rewritten as 

i J d^kduje^^^e-^^^-^iuj -k-V- I'ujc - ivNL)5fv = Sr o J d^kduje^^^e-^^^'^5^U ■ V/o -f 

Equation (I149p is always the Boltzmann equation applied to 6f and described in the Fourier dual space. The 
gyrokinetic equation is the same equation with 1 ' = 0 : 

i J d^kduje^‘^*'e~^^'^ {(jj — k-V — WML)dfo = 

= dtSfo + V-VSfo + iyNLSfo = J d^kduje^‘^^e-^'^-^S^U-Vfo + 


(150) 


-i / d^kduje^^*e-^'^-^ — [m-V-SA)Jo + ^^bxSA-kJi]dJo, 
J m \kt\ 


where 


Sfo = 




(151) 


is the gyro-average of Sf and it represents the gyro-centers in Sf. The solution for Sfp is 

. ~ [(<5$ - E • dA)Ji, - {pLu;c/\kt\)k x b ■ 6Aj;,]dJo 

ojl' =- ] - 77 - 7, -^- 

m UJ — k • V — I'UJc — Ik'NL 

t6i,oe^>^<S^U-Vfo 

oj — k ■ V — I'ujc — ivNL 

In particular, for V = 0, the solution 5/o of the gyrokinetic equation is 

[( 5 $ _ y . 5A) Jo + {pLOOc/\kt\)k x b ■ SAji]d,fo , 


(152) 


Sfo = -(ew/m)- 


cj — k ’ V — ivnl 


—i 


e^’^<5^U ■ V/o 
OJ — k ■ V — iuNL 


(153) 













20 


In real space, the solution of the Boltzmann equation is 

OO « 

5f = -Y^ / ( 154 ) 

1 — — 00 

[((5$ - V ■ SA)Ji - {pLUJc/\kt\)k X b ■ SAJl]defo + i{m/e)5iQe^^'^{5(_U/w) ■ V/p 

Lo — k ■ V — IlOc — ivNL 

Such solution is formal in the sense that it employes the symbol vnl which is not a well defined mathematical object. 
If C is known, then the problem could be attacked with the (perturbative) functional analysis and, maybe, solved. 
There is another limiting behavior of the system in which the solution can be given without explicitly knowing the 
collisional operator. This happens at low collisionality (if (7 = 0 the equation (11331) is said collisionless Boltzmann 
equation) when the magnetic field is slowly varying and the initial distribution function /o is an equilibrium one: 
fo feq- In such case the system is said to be at marginal stability. In this case, the distribution function Sf is very 
small if compared to feq] moreover, F —>■ C/, or ^ 0, and —>■ iz becomes a real number. At marginal stability, 

the solution ()153|) for the gyrokinetic equation (I150p is well defined (if n replaces and ^ = 0): 


Sfo \m.s.= -(eujlm) 


[((5$ - U ■ 5A)Jo + {pLaJc/\kt\)k x b ■ SAji]dsfe, 


uj — k ■ U — iv 

The same is true for 1)15411 which is solution of the Boltzmann equation at marginal stability: 

[(5$ - U ■ SA)Ji> - {pLUJc/\kt\)k x b ■ SAJl,]defeq 


5fi' 


m 


(jj — k ’ U — VuJc — iiy 


(155) 


(156) 


In the general case, the differences between p53L where the quantities are not obtained perturbatively, and the 
standard gyrokinetic equation, e.g. 


22 ) 1 . are mildly overcome if defo is written explicitly: 


defo = deSo ■ dejo + d^Q ■ V/o + dePo ' df,Jo, (157) 

and, after that, only the first order quantities are considered: b —>■ b^^\ U uy and ojc —t (e|i?|/m). However, a one 
to one correspondence seems to be lost and it should be further investigated. Indeed, such differences can be very 
important for stability studies of plasmas, as recently reviewed in (^ . 


X. CONCLUSIONS 

In the present work it has been shown the possibility of analytically describing the non relativistic motion of charged 
particles in general magnetic fields. The method used to address this problem is the same approach used by the gyro- 
center transformation but adopting a non perturbative magnetic field description to find new solutions of motion in 
GC or gyro-center coordinates. The equation (HU has been easily obtained from combining the main equation for a 
dynamic system, o, the Lorentz force law and the condition (USD. From the equation (Ha two solutions have been 
proposed: the GP and the gyro-particle solutions. The first has been characterized to be minimally coupled to the 
magnetic field, the second has been defined to be maximally coupled to the magnetic field and also to move on a closed 
orbit. A generic charged particle has been shown to have a velocity expressed as the sum of such particular solutions, 
dZZl) and dUl), and characterized by the gyro-centers coordinates {X,V), a Larmor radius pL and a gyro-angle 7 . 
The GC description of motion is similar to the gyro-center description but with the gyro-center coordinates A, V 
computed at a particular time, here settled at t = 0. The explicit motion description has been shown in section EH 
where the magnetic field has been chosen to have an axial symmetry. The axisymmetric toroidal magnetic field, often 
encountered for describing laboratory plasmas, has been extensively studied obtaining new relevant representations 
for the velocity of particles and for the toroidal canonical momentum. The symmetry in the gyro-angle has revealed 
the importance on having defined the gyro-particle as moving on closed orbit. Such condition ensures the existence of 
a particular gauge function (flMll . which is proportional to the product of the magnetic moment of the gyro-particle 
times the gyro-angle. In such gauge, with the defined coordinates, it has been possible to represent the magnetic 
field seen by the particle as an axisymmetric magnetic field. The gyro-angle 7 becomes an ignorable coordinate and 
the Larmor radius is shown to be proportional to the magnetic flux linked to the closed orbit. It is furnished the 
definition (ITT71) of the magnetic moment in gyro-center as in GG coordinates. The magnetic moment is considered an 
exact constant of motion instead of an adiabatic one, and also the drift velocity is expressed in an exact form (11251) . 
Such drift velocity representation is particularly simple because it is very similar to the drift velocity developed at the 
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first order in the standard gyro-center transformation. Indeed, it is shown to be constituted by the E x B-like drift, 
the V|i?|-like drift, the curvature-like drift, the i9t|i?|-like drift and the i?anos-like drift. Once a particular set of GC 
coordinates (t, X, e, /r, 7 ), is chosen to describe the motion of a charge in a general magnetic field, the first use of the 
present approach is on deriving the gyrokinetic equation. Such derivation of gyrokinetic equation from the Boltzmann 
equation without explicitly considering the collision operator, is shown without any approximations (neither ordering 
considerations). The final form (115011 . which is particularly simple, is formally solved, (I155L together with the formal 
solution, (ITMl) . of the Boltzmann equation for the perturbed distribution function Sf. The obtained solution can 
be considered the correct solution for the Boltzmann equation at marginal stability. The next steps, concerning the 
use of the present non perturbative approach to gyrokinetics, should be to address the Vlasov equation, where the 
fluctuation of fields is coupled to Sf via Maxwell equations, and to study its phenomenologies, e.g. obtaining the 
plasma dispersion relation. 
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